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Abstract 

A Skyrme-type momentum- dependent nucleon-nucleon force distinguishing isospin effect is 
parameterized and further implemented in the Lanzhou Quantum Molecular Dynamics (LQMD) 
model for the first time, which leads to a splitting of nucleon effective mass in nuclear matter. 
Based on the isospin- and momentum- dependent transport model, we investigate the influence 
of momentum-dependent symmetry potential on several isospin-sensitive observables in heavy-ion 
collisions. It is found that symmetry potentials with and without the momentum dependence but 
corresponding to the same density dependence of the symmetry energy result in different distri- 
butions of the observables. The mid-rapidity neutron/proton ratios at high transverse momenta 
and the excitation functions of the total 7r~//T + and K°/K + yields are particularly sensitive to 
the momentum dependence of the symmetry potential. 
PACS: 24.10.Lx, 25.75.-q, 25.70.Pq 
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Heavy-ion collisions induced by neutron-rich nuclei at intermediate and relativistic energies in 
terrestrial laboratories are a useful tool to extract the information of nuclear equation of state 
(EoS) of isospin asymmetric nuclear matter, which is poorly known for the high-density neutron- 
rich matter but has an important application in astrophysics, such as the structure of neutron 
star, the cooling of protoneutron stars, the nucleosynthesis during supernova explosion of massive 
stars etc [TJ. The EoS of nuclear matter is usually expressed through the energy per nucleon as 
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E(p, 5) = E(p, 5 = 0) + E S y m (p)5 2 + 0(S 2 ) in terms of baryon density p = p n + P P , relative neutron 
excess 5 = (p n — p p )/(p n + P P ), energy per nucleon in a symmetric nuclear matter E(p,5 = 0) 
and bulk nuclear symmetry energy -Esym = \ d f^' 5 ' 1 1 5=0- Based on several complementary 
analysis of available experimental data associated with transport models, a symmetry energy of 
E S y m (p) ps 31.6(p/p ) 7 MeV with 7 = 0.69 — 1.05 was extracted for densities between 0.1p and 
1.2po |2]. However, predictions for high-density symmetry energies based on various microscopical 
or phenomenological many-body theories diverge widely [3j HI [5]. More realistic approach to 
extract the information of the -Esym(p) is the comparison between transport model calculations 
and experimental data. For that, the reliable input potentials including the isovector (symmetry 
potential) and isoscalar parts in the transport models are very necessary for precisely predicting 
some promising observables. 

The momentum dependence of the isoscalar potential leads to the same nucleon effective 
mass (Landau mass) for neutrons and protons in nuclear matter and has been widely studied in 
transport models for heavy-ion collisions. The effective mass splitting of neutrons and protons 
results from the momentum-dependent symmetry potential, which has been implemented in the 
one-body transport models, such as the isospin Boltzmann-Uehling-Uhlenbeck (IBUU04) model 
[6] and the stochastic mean-field (SMF) model In this work, a Skyrme-type momentum 

dependent nucleon-nucleon (NN) force distinguishing protons and neutrons is parameterized and 
then included in a N-body approach (LQMD model) for the first time. Furthermore, effects of the 
momentum dependence in heavy-ion collisions are investigated. In particular, its influence on the 
isospin sensitive observables to extract the high-density symmetry energy is discussed. 

The LQMD model has been successfully applied to treat the dynamics in heavy-ion fusion 
reactions near Coulomb barrier and also to describe the capture of two heavy colliding nuclides to 
form a superheavy nucleus [9J, [10] . Further improvements of the LQMD model were performed in 
order to investigate the dynamics of pion and strangeness productions in heavy-ion collisions and 
also to extract the information of isospin asymmetric EoS at supra-saturation densities [TTJ [TSJ 
[131 [H]. I n the previous versions, we only considered the scalar part of the momentum-dependent 
interaction and the density-dependent symmetry potential. We have included the resonances 
(A(1232), N*(1440), N*(1535)), hyperons (A, E) and mesons (n, K, 77) in hadron-hadron collisions 
and the decays of resonances for treating heavy-ion collisions in the region of 1A GeV energies. 
The time evolutions of the baryons (nucleons and resonances) and mesons in the system under 
the self-consistently generated mean-field are governed by Hamilton's equations of motion, which 
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read as 

OH OH 

We only consider the Coulomb interaction for charged hyperons. The Hamiltonian of baryons 
consists of the relativistic energy, the effective interaction potential and the momentum dependent 
part as follows: 



R B=Y. "JPi + m i + U int + U m om. (2) 
% 

Here the Pj and rrii represent the momentum and the mass of the baryons. 

The effective interaction potential is composed of the Coulomb interaction and the local inter- 
action 

Uint = Ucoul + U[ oc . (3) 

The Coulomb interaction potential is written as 

Ucoui = ~ E ^erffa/VU) (4) 

where the €j is the charged number including protons and charged resonances. The = |r$ — r^| 
is the relative distance of two charged particles. 

The local interaction potential is derived directly from the Skyrme energy-density functional 
and expressed as 

Ui oc = J V loc {p{r))dv. (5) 
The local potential energy-density functional reads 

vUp) = f P - + t^~^~t + f^(Vp) 2 + |^[v(p„ - Pp )f 

2 po 1 + 7 Po 2p 2p 

+E% c m {p)p5 2 + g T p^/pl l \ (6) 

where the p n , p p and p = p n + p p are the neutron, proton and total densities, respectively, and 
the 5 = (p n — Pp)/(p n + Pp) is the isospin asymmetry. The coefficients a, (3, 7, g sur , gf° r , g T are 
related to the Skyrme parameters to, £1,^2, £3 and xo, x±, £2, £3 [ID] . The E l s ° c m is the local part of 
the symmetry energy, which can be adjusted to mimic predictions calculated by microscopical or 
phenomenological many-body theories and has two-type forms as follows: 

E% c m { P ) = \c sym {p/ Po y% (7) 



and 



El s°vm(p) = a Sym (p/po) + b sym (p/pof 



sym\ 



The parameters C sym , a sym and b sym are taken as 52.5 MeV, 43 MeV, -16.75 MeV and 38 MeV, 37.7 
MeV, -18.7 MeV for the cases with and without momentum- dependent interactions, respectively. 
The values of 7 S =0.5, 1., 2. correspond to the soft, linear and hard symmetry energy, respectively, 
and the Eq. (8) gives a supersoft symmetry energy, which cover the largely uncertain of nuclear 
symmetry energy, particularly at the supra-saturation densities. 

We have taken the same Skyrme-type form for the momentum-dependent potential in the 
Hamiltonian as in Ref. [15] but distinguishing isospin effect, which is expressed as 



Umom = ^— J2 C r,r' 5 T,rA',T 3 J J J dpdp dr f { (r, p, t) [ln(e(p - p') 2 + 1)] 2 fj (r, p', t) . (9) 

™ t,t' 

The term is also given from the energy- density functional in nuclear matter, 

U mom = ^-J2 C ry J J /rfprfp , rfr^(r, P )[ln( e ( P - P ') 2 + l)] 2 ^'(r,p , )• (10) 

' " t,t' 

Here C T)T = C mom {l + x), C T y = C mom {l — x) (t / r') and the isospin symbols r(r') represent 
proton or neutron. The sign of x determines different mass splitting of proton and neutron in 
nuclear medium, e.g. positive signs corresponding to the case of m* < m*. The parameters 
Cmom arid e were determined by fitting the real part of optical potential as a function of incident 
energy from the proton-nucleus elastic scattering data, which determine the nucleon effective 
mass in isospin symmetric nuclear matter. In the calculation, we take the values of 1.76 MeV, 
500 c 2 /GeV 2 and -0.65 for the C mom , e and x, respectively, which result in the effective mass 
m* /m—0.75 in nuclear medium at saturation density for symmetric nuclear matter. For the cold 
nuclear matter, we have the phase-space density / r (r, p) = p T (r)<d(p F (r) — |p|)/(47T^(r)/3) with 
the Fermi momentum Pf( t ) — h(3ir 2 Pt) 1 / 3 . Implementing the phase-space distribution into Eq. 
(10), we get the contribution of the momentum dependence to the symmetry energy E™^(p). 
Therefore, the symmetry energy per nucleon in the LQMD model is composed of three parts, 
namely the kinetic energy, the local part and the momentum dependence of the potential energy 
as 

1 h 2 /3 \ 2 / 3 

E sym (p) = 3^ y?P) + E% c m{p) + ET y T(p)- (ii) 
Figure 1 is a comparison of different stiffness of nuclear symmetry energy after inclusion of the 
momentum-dependent interactions. One can see that all cases cross at saturation density with 
the value of 31.5 MeV. The local part of the symmetry energy can be adjusted to reflect the 
uncertain behavior of the symmetry energy at sub- and supra-normal densities. The contributions 
of the local and momentum-dependent interactions of the potential part in the symmetry energy is 
shown in Fig. 2. The momentum dependence has a positive contribution to the symmetry energy 
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in the mass splitting of m* n < m*, which is opposite for the case of m* > m*. The local part of 
the symmetry energy is adjusted to get the same density dependence of the symmetry energy for 
a given stiffness with increasing density. 

Combined Eq. (6) and Eq. (10), we get a density, isospin and momentum- dependent single- 
particle potential in nuclear matter as follows: 

u T ( P , 5, P ) = a £. + fj£ + \ 9tP ^/pT + E l ° c m (pW + ^^ P 6 2 + E\- m { P ) P f- 

Po Po <j up op T 

+ IC [ rf P 7 T (r,p)[ln( e (p-p') 2 + l)] 2 
Po J 

+ -C T y f dp'^(r,p)[ln(e(p - p') 2 + l)] 2 . (12) 

Here r ^ t', d5 2 /dp n = 45p p /p 2 and d5 2 /dp p = —45p n /p 2 . The effective (Landau) mass in nuclear 
matter is calculated through the potential as m* = m T / (l + jpjl^p |) with the free mass m T at 
Fermi momentum p = p F . Therefore, the nucleon effective mass only depends on the momentum- 
dependent interactions. The mass splitting of protons and neutrons in nuclear matter as functions 
of density (5=0.2) and isospin asymmetry (5 = (p n — p p ) / \p n + p p )) at saturation density is shown 
in Fig. 3. The left windows are the cases of m* < m* with the parameter x=0.65 in the coefficient 
C T y and the right panels are m* > m* The two sorts of the mass splitting can be chosen in the 
LQMD calculations. In accordance with the Lane potential [IB], the symmetry potential can be 
evaluated from the single-nucleon potential U sym (p, p) = (U n (p, S, p) — U p (p, 5, p))/2S. 

A hard core scattering in two-particle collisions is assumed in the simulation of the collision 
processes by Monte Carlo procedures, in which the scattering of two particles is determined by a 



geometrical minimum distance criterion d < J 0.1 <Jtot/^ fm weighted by the Pauli blocking of the 
final states [T7JHE]- Here, the total cross section a tot in mb is the sum of the elastic and all inelastic 
cross sections. The probability reaching a channel in a collision is calculated by its contribution of 
the channel cross section to the total cross section as P c h = cy c h/a tot . The choice of the channel is 
done randomly by the weight of the probability. The primary products in nucleon-nucleon (NN) 
collisions in the region of 1A GeV energies are the resonances A(1232), iV*(1440), iV*(1535) and 
the pions. We have included the reaction channels as follows: 

NN <B- iVA, NN <-> NN*, NN <B- AA, 

A ^ A%, N* <r> Nti, NN -> NN7i(s - state), iV*(1535) -> Nr]. (13) 



At the considered energies, there are mostly A resonances which disintegrate into a it and a 
nucleon in the evolutions. However, the iV* yet gives considerable contribution to the energetic 



pion yields. The energy and momentum-dependent decay widths are used in the calculation [TTJ 
for the A(1232) and A r *(1440) resonances. We have taken a constant width T=150 MeV for the 
iV*(1535) decay. The strangeness is created by inelastic hadron-hadron collisions [T4J [19]. We 
included the channels as follows: 

BB ->• BYK, BB ->• BBKK, Bn -> YK, Bn -> JV/nT, 

Yn ->■ £if , -> Ftt, FiV -> ffl (14) 

Here the B strands for (N, A, N*) and Y(A, £), K(K°, K+) and F(7f°, K - ). The elastic scattering 
between strangeness and baryons are considered through the channels KB — » fTi?, K.B — >■ YB 
and — >■ A'i?. The evolutions of mesons (n, K, 77) are also determined by the Hamiltonian in 
which the Coulomb interaction and the in-medium potential were considered in the model [131 EE] • 

To check the influence of the momentum dependence of the symmetry potential on reaction 
dynamics, we calculated the transverse momentum distributions of the ratios of neutrons over 
protons in the mid-rapidity domain and the excitation functions of charged pion yields as shown 
in Fig. 4. Inclusion of the momentum-dependent interaction in the symmetry potential reduces 
the n/p yields at high transverse momentum owing to its negative contribution to the symmetry 
energy at the case of m* > m* in nuclear medium, which enforces an attractive force in neutron- 
neutron collisions and further increases the collision probabilities. Contrarily, the case without the 
isovector of the momentum-dependent interaction quickly squeezes out neutrons, in particular for 
the energic neutrons, which nearly appears a flat distribution at p t >0.3 GeV/c. The conclusions 
are consistent with the calculations based on the IBUU04 transport model [6]. The charged pion 
ratios are slightly changed by the momentum-dependent potential. Therefore, the neutron and 
proton transverse emission ratio in the mid-rapidity region is to be a nice probe of the isovector 
part of the momentum-dependent interaction. 

The pions and strange particles are mainly produced in the dense hadronic matter formed 
in heavy-ion collisions. The tt~/tt + and K°/K + ratios can be probes to extract the high-density 
information of isospin asymmetric EoS. Shown in Fig. 5 is a comparison of the excitation functions 
of the 7r~/7r + and K°/K + ratios with and without the momentum dependence of the symmetry 
potential in the 197 Au+ 197 Au reaction for head-on collisions. One can see that the momentum- 
dependent interaction by distinguishing isospin effect reduces the 7r _ /tt + and K°/K + yields, which 
is more pronounced close to the threshold energies of meson production. The results are caused 
from the fact that the interaction enhances the energic neutron-neutron collisions. Furthermore, 
the 7T~ and K° is produced through the channels nn — > pA~ , A~ — > n-n^ and nn — > nAK°. 



In summary, influence of the momentum dependence of the symmetry potential on isospin 
sensitive observables in heavy-ion collisions is investigated by using an isospin- and momentum- 
dependent transport model (LQMD). It is found that the momentum dependence of the symme- 
try potential plays an important role on nucleon transverse emissions and the ratios of 7i^/n + 
and K°/K + , which are also as promising probes of high-density symmetry energy. To precisely 
constrain the density dependence of the nuclear symmetry energy, one firstly needs to get the 
accurate information of the density- and momentum-dependent symmetry potential. Further ex- 
perimental data of determining the momentum-dependent symmetry potential are very necessary 
to investigate accurately the dense neutron-rich matter. The updated LQMD model would to be 
a useful tool to predict the density dependence of the symmetry energy from heavy-ion collisions, 
in particular at supra-normal densities, and also to analyze experimental data for constraining the 
symmetry energy. 
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Figure 1: Density dependence of the nuclear symmetry energy for the cases of supersoft, soft, 
linear and hard trends. 



Q 



asym-hard asym-supersoft 




p/p p/p 



Figure 2: Density dependence of the potential part of nuclear symmetry energy at different mass 
splitting for hard (left panel) and supersoft (right panel) symmetry energies. 
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Figure 3: Nucleon effective mass normalized by vacuum mass as functions of density (5 = 0.2) 
and isospin asymmetry (p = p ). 




p t (GeV/c) p t (GeV/c) 



Figure 4: Transverse momentum distributions of the ratios of neutron/proton within the rapidity 
bin \y/y pro j\ <0.3 and n~ /n + in central 124 Sn+ 124 Sn collisions at incident energy 400A MeV. 
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Figure 5: Comparison of excitation functions of the n /ir + and K° / K + yields for head-on collisions 
in the 197 Au+ 197 Au reaction with and without the isovector part of momentum dependence. 
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